*------------------------------------------------------------------------------*
				*** FIGURE S10. BOXPLOT PM10 WEEK vs WEEKEND ***
*------------------------------------------------------------------------------*

use PM10_ug day month year using "${tempdata}pollution_UBA_d_19902019_w_av.dta", clear

keep if year >= 2000 & year <= 2018

*------------------------------------------------------------------------------*

*Divide between week days and weekend
gen dow = dow( mdy( month, day, year) )
gen PM10_10ug = PM10_ug / 10
gen weekend = 0
replace weekend = 1 if dow == 6 | dow == 0

*------------------------------------------------------------------------------*

***Box-Plot
graph box PM10_10ug, over(weekend, relabel(1 "Monday-Friday" 2 "Weekend")) nooutsides note("") ytitle("PM10 (10 {&mu}g/m{sup:3})") b1title("Days of the week") graphregion(color(white))

graph export "${outputs}figure_S10.pdf", replace

*------------------------------------------------------------------------------*

clear

exit
